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ABSTRACT 

By redshift of 10, star formation in the first objects should have produced 
considerable amounts of Carbon, Nitrogen and Oxygen. The submillimeter lines 
of C, N and O redshift into the millimeter and centimeter bands (0.5 mm-1.2 cm), 
where they may be detectable. High spectral resolution observations could potentially 
detect inhomogeneities in C, N and O emission, and see the first objects forming 
at high redshift. We calculate expected intensity fluctuations and discuss frequency 
and angular resolution required to detect them. For CII emission, we estimate the 
intensity using two independent methods: the line emission coefficient argument and 
the luminosity density argument. We find they are in good agreement. At 1 + z ~ 10, 
the typical protogalaxy has a velocity dispersion of 30kms~^ and angular size of 
1 arcsecond. If CII is the dominant coolant, then we estimate a characteristic line 
strength of ~ O.lKkms^^. We also discuss other atomic lines and estimate their 
signal. Observations with angular resolution of 10"'^ can detect moderately nonlinear 
fluctuations of amplitude 2 • 10~^ times the microwave background. If the intensity 
fluctuations are detected, they will probe matter density inhomogeneity, chemical 
evolution and ionization history at high redshifts. 
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1. Introduction 



When did tiie first stars and galaxies form? How can we detect them? Quasars were once 
the most distant objects that have ever been observed ( Schneider, Schmidt and Gunn 199l| ). Now 
galaxies are beginning to take their place. Discovery of a large number of galaxies at z > 3 ( ^teidel 
et al. 199^ ), and even up to z ~ 5 ( Franx et al. 1997 , Dey et al. 1998| ), suggests that galaxies existed 
well before quasars did. In hierarchical models, star formation starts relatively early, by z ~ 30 



( Couchman and Rees 198t , Fukugita and Kawasaki 1994 , Gnedin and Ostriker 199^ , Tegmark et 
al. 1997), so even these z ~ 5 galaxies are not representative of the first generation of objects. 



The purpose of this paper is to investigate the possibility of detecting primordial objects 
through the atomic lines of carbon, nitrogen and oxygen, which should have been produced as a 
result of star formation. At low redshift, these submillimeter lines are important coolants (e. g.. 



Bennett et al. 1994 




Israel, White and Baas 1995 


et al. 1997, Malhotra et al. 1997). If the source is 



redshift 1 + z ~ 10-20, the lines are redshifted 
to radio frequencies, where the atmosphere is more transparent and ground-based detectors with 
very high sensitivity and resolution exist. 

If the source was uniformly distributed through space, then its emission would merely distort 
the spectrum. Since heavy elements produced by stars are distributed inhomogeneously, the 
observed emission should have intensity fluctuations as a function of both position and frequency. 
With high enough resolution one can in principle distinguish them from background. If we take 
Ih-^ Mpc characteristic comoving scale of the inhomogeneity, this corresponds to ~ 1 arcmin 
in angular scale, and ~ 10^'^ in frequency resolution. On smaller scales, velocity dispersion of 
matter limits the scale that can be resolved in frequency space. In this paper we formulate the 
intensity fluctuations and estimate their magnitude. This signal should give us information on 
number density fluctuations of specific source elements at 1 + z > 10. Detailed determination of 
the peak of this signal will also tell us at which epoch elements were ionized. 

The rest of this paper is organized as follows. In the next section, we present our formalism 
for estimating the fluctuating signal. Then we estimate its magnitude in the case that the matter 
density fluctuations are calculated from the linear power spectrum in section |3|. We also estimate 
optimal frequency and angular resolution for an object search. We focus on atomic line redshifted 
into present radio band, investigating carbon, nitrogen and oxygen emission lines at rest frame 
wavelength of hundreds of micrometers. In section ^, we discuss the size and the number density 
of nonlinear clumps at the reionization epoch, and estimate intensity from nonlinear clumps and 
the interval between enhance of intensity fluctuations. Section |5| gives discussion of merits and 
demerits as the signal of specific emission. After alternative estimation of CII luminosity in section 
^, the final section gives our conclusion. 



2. Observable intensity fluctuations 
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2.1. Formulation 

Throughout this paper we make the following assumptions: by 1 + z ~ 10-20, star formation 
has both formed metals and heated the surrounding gas initially to T ~ 10^ K. Cooling is 
compensated by heating so that the gas temperature Tgas is maintained much higher than the 
CMB temperature at that epoch but lower than ~ 10^ K. Cooling by atomic transitions exceeds 
Compton cooling, as we will see below. 

Intensity luiz) of radiation in emitting medium at redshift z is obtained by solving the 
radiative transfer equation: 

^^ = ci,(z)-3F(z)J,(z), (1) 

where c is the speed of light, and juiz) and H{z) denote the spontaneous emission coefficient and 
Hubble parameter, respectively. Differentiation with respect to time in the left hand side is taken 
with + z)~^ = const. The second term in the right hand side comes from the effect of cosmic 
expansion. We have assumed that absorption is negligible. Integrating equation (||) yields the 
present intensity at frequency i^q: 

where ly = + z) in the integrand. For spontaneous transition from an upper level n to a lower 
level I releasing a photon with energy /ifunej the emission coefficient is 

Ju = nuAui(p{iy), (3) 

where hp is the Planck constant, Uu is the number density of the source in the upper level, Aui 
is the Einstein A-coefficient, and (piiy) = Soii^ — z^Une) is the line profile function, which we take 
to be equal to delta function because both natural and thermal broadening is negligible for the 
expected frequency resolution. Since redshift of the source Zg is inferred from the relation: 

l + ^s = ^, (4) 

we can express the line profile function as 

(t>{iy) = (t>{M^ + z)) = ^^6D{z-z,). (5) 

Using equations (^) and (^) in equation (|2|) yields intensity due to the line emission: 

hne = ^AuinuH-\z,){l + z,)-\ (6) 
47r 

In the Einstein-de Sitter universe, equation (P) reads 

hne = ^H^'A^iUn,{l + z,y^/\ (7) 
47r 

where = riu/ns^ is the fraction of the source that is in the upper state, and Ua is the local 
number density of the source. 
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2.2. Magnetic dipole transition 

The spin-orbit interaction splits each spectroscopic term of an atom into a set of energy 
levels. Magnetic dipole transitions near the ground state of neutral and ionized carbon, nitrogen 
and oxygen have wavelength of hundreds of micrometers. We list such transitions in table || 
(Melnick 1988). These are candidates which produce intensity fluctuations. 



2.3. Critical number density 

We have seen the line intensity is proportional to the number density of the source that is in 
the upper state. Let us examine the condition for substantial number of the source to be in the 
upper state. If either spontaneous or stimulated emission occurs more rapidly than the collisional 
excitation by surrounding gas (mainly Hydrogen), very few atoms (or ions) will be in the upper 
state. The condition for collisional rate exceeding spontaneous and stimulated emission rates is 

ngas > (8) 
where rig^^g is the critical number density determined by 



1 — exp 



kTcMB{Zs 



(9) 



with 7„/ being the rate coefficient for collisional de-excitation, and Tcmb(-z) the CMB temperature 
at redshift z. 

If Tigas > "-gls) the excitation temperature of the source Tgource is almost equal to Tgas, where 
Tgas is maintained around 10^ K. In this case, the fraction fu = does not depend on Tgource- 
In CI case, the excitation energies of and levels correspond to 60K and 1.5 x lO^K, 
respectively. If we assume that the gas is in thermal equilibrium at a temperature around lO^K, 
electrons in the (2p)^ shell of CI easily transit within levels, while it is almost impossible to 
jump to level. Therefore the fraction of atom that is in the state labeled as u is independent 
of matter temperature and determined by 

= ^ = , (10) 

'^tot Z^same term 9i 

where = 2 J + 1 is the statistical weight of level u. Specifically, the fraction is 1/9, 3/9 and 5/9 
for ^Pq,^Pi and of neutral carbon. 

If f^gas < "-gis ) the line intensity is diminished to 



(11) 
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In CI (609 iJ,m) case, the critical number density of Hydrogen is given by 



73 



103 K 



-0.34 



1 



exp 



8.7 

1 + Zs 



-1 



cm 



(12) 



(Hollenbach and McKee 1989). The last factor in the above expression comes from the effect 
of stimulated emission. The effect is important for the CI 609 fim line, but is negligible for the 
other fine structure transitions of interest here at 1 + ^ 20. The excitation of ions is induced 
mainly by the collisions with electrons, if the gas is ionized. The condition (^) is determined by 
the critical number density of electron 



28 



103 K 



0.50 



cm 



(13) 



([Hollenbach and McKee 1989|) , in CII (158/im) case. 



2.4. Gas clumping 



At 1 + z ~ 10-20, the mean density of gas 



2 X 10 [(1 + z)/10]3 cm 3 is lower than 



the critical number density. The condition for sufficient magnitude of signal, however, is attached 
to local number density of the source. The condition (@) can hold in the virialized object, even if 



Let us express by F'^^ the fraction of gas within the relevant scale which satisfies Ug 
Then the smoothed intensity is given by 



> n 



June = ^H^'Auirj^F^ 



nH,o(l + 2s) 



-3/2 



gas- 



(14) 



where [na/rin] is relative abundance of the atom, and nH,o is present Hydrogen number density. 
It is uncertain to what degree the baryonic gas condenses in dark matter halos. However, it 
seems reasonable to expect that the gas inside such halos can radiatively cool and form a gaseous 
disk, whose radius is Arvir, where A is the spin parameter and rvir is the virial radius ( Peebles 



1993, Dalcanton, Spergel and Summers 1997). If this is the case, we expect the mean gas density 
to be at least 187r^A~3 (~ 10^) times the mean cosmic gas density at collapse, where A = 0.05 



(Barnes and Efstathiou 1987) typically. For objects that form at 1 -|- z > 10, this implies a 
characteristic number density of baryonic gas of at least 300cm~3. With this consideration in 
mind, is estimated by the fraction of matter within the collapsed objects with mass M greater 
than the Jeans mass Mj. By 1 -|- z = 10, the fraction amounts to 50 percent of the total mass 
in the universe ( Fukugita and Kawasaki 1994 ), assuming Mj ~ IO^Mq. In the case that gas is 
considerably reheated, Mj rises to ~ IO^Mq and = 0.3 ([Fukugita and Kawasaki 1994]) . 



Gas clumps in these halos are optically thin in the relevant lines, and therefore the formulation 



in Section 2T is valid. For example, the one-dimensional velocity dispersion of a W^Mq halo 
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at 1 + z = 10 is estimated to be ~ 3kms~^ using the spherical cohapse model. The size of 
the gas clump inside the halo with rigas = 10^ cm~^ is ~ lO/i^^pc. With these numbers and 
'^c/^H = O.Ol(n(;;;/'raH)0, the optical depth in the CI 609 /zm line is estimated to be around 10~ 



Intensity fluctuations with frequency correspond to number density fluctuations around 
relevant redshift: . , « , 

^= ^) . (15) 

-'line \ ^ / source 

Even if the intensity /une itself is smaller than the CMB intensity IcmBj the fluctuations 
A/iinc/-^tot — A-^iinc/-^CMB Can be extracted as difference between those at different frequencies. 



3. Intensity fluctuations due to linear perturbations 

In this section, let us estimate the signal originated from the inhomogeneities on a 
characteristic scale ~ lh~^ Mpc corresponding to linear perturbation. In the next section, we 
focus on smaller scales and estimate the intensity from nonlinear clumps. 



3.1. Magnitude of the signal from carbon 609 //m emission 

In this and next sections, we take CI (609 ^m) as an example ( [Yamamoto and Saito 199 Ij ). 
Essential discussion is common with other lines. Since the ionization state of the universe at z > 5 
is not known, we begin by assuming that most of the carbon at 1 + z = 10 is in the form of CI. 
The ratio of the intensity of CI (609 fxm) line emitted at 1 + 2 = 10 to total intensity /tot can be 
calculated from equation (14) with f^^ = 1/3 (Section |2.3|) : 



In, 



'tot 



ICMB 



(16) 



V0.5/ \^0.025y VO.Ol(nc/nH)0/ V3.6 X 10-V Vo.Sy ' ^ ^ 

where h is the present Hubble parameter in units of 100 kms^^Mpc^^, and JIb is the present 
baryon density parameter. It is reasonable to expect that nc/nu = O.Ol(nc/nH)0 by 1 + z = 10. 
Numerical simulation ( Pnedin and Ostriker 199^ ) of nonlinear self-gravitating clumps at high 
redshift shows they have a metallicity Z = O.O3Z0 in the high density regions. Observation of CIV 
absorption in Ly a forest clouds also supports the above value of nc/nu. In this paper, we adopt 
(nc,0, reN,0, J^o,©) = (3.58, 1.12,8.49) x lO^^nn,© for the solar abundance of carbon, nitrogen and 
oxygen ([Anders and Grewesse 1989| , [Woosley and Weaver 1995 ). 



For 1 + z ^ 15, if most of the Carbon is in its ground states and ngas > ?T-gis, then cooling by 
the CI line emission exceeds Compton cooling: the cooling time due to the CI line emission can 



-7- 



be estimated as 

where /cb is Boltzmann's constant. On the other hand, Compton coohng time is ( |Levin, Freese| 
and Spergel"T99^ ) 

yr, (19) 



^8 



^Compton ~ 1 X 10' 

where rig is the number density of free electrons. 



, rip 
1+ ( — 



10 



We can estimate number density fluctuations in the standard cold dark matter (CDM) model 
with h = 0.5, the density parameter JIq = Ij and the present density fluctuations on 8/i~^ Mpc 
scale cjg = 1: 

(An/n),tom = AM/M = 0.4 (20) 

on l/i^^ Mpc scale at 1 + z = 10. This has been calculated from the CDM power spectrum P{k) 
given by Bardeen et al. (1986). This is a conservative estimate as the predicted fluctuation level is 
even larger in vacuum-dominated or open models. Therefore the magnitude of the signal becomes 
A^iinc/-^CMB = 1 X 10^^ in the standard CDM model. Table Q summarizes predictions of intensity 
fluctuations with other elements. Note that each value has been calculated assuming that 100% 
of the element is in the form of the neutral atom or the ion in the first column. Actual signal is 
reduced by the fraction of the species in each element. 



Since the metals are produced in dense regions, they will be "biased" tracers of the underlying 
mass. This implies that 

\ / atom 

where h is the bias parameter. Gnedin and Ostriker (1996) find a strong correlation between 
metallicity and density. If Z oc p", then in the linear regime, this implies h = 1 + a, which enhances 
the signal over the values estimated in table |l[ 



3.2. Characteristic fluctuation size 

Now, we investigate the resolution that we need to detect the intensity fluctuations due to the 
line emission. The comoving length scale Rs probed with frequency resolution Ai//i/ is given by 

A c(l + z) Aa c(l + z) Az/ 

Rs = cAr] = = , (22) 

H[z) a H{z) V 

where r\ is the conformal time, and a = (1 + z)^^ . In the Einstein-de Sitter universe, 

R, = cH^\l + z)-^I^Au/u = 0.95 h-' Mpc. (23) 
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An observation with a resolution jv means that we probe the matter density inhomogeneities 
smoothed by spatial thickness i?s- Figure |l| shows predicted smoothed intensity fluctuations in 
the standard CDM model which are observed with given resolution Ai//i/ in CI (609 /um) and CI 
(370 /xm) cases. It also shows the signal for CII (158 ^m) in the case where most carbon is in 
the form of CII. Higher frequency resolution brings larger signal because of larger matter density 
fluctuations. For appropriate correspondence of frequency space to source positions, iig is required 
to satisfy the relation: 

^(i?3) > (24) 
V c 

Here i'random(-Rs; -2) is the random velocity inside a spherical region of comoving radius R^,: 

^^random(fis; z) = J dkP{k) [l - W^R^; k)\ , (25) 



where W{Rs; k) = 3(sin kR^ — kRg cos kRg) /{kR^}^. In the standard CDM model, the relation ( p4D 
yields Au/u > 2 x 10"^ and i?s > 0.02/1"^ Mpc. 

Let us move on to the needed angular resolution. To maximize observed intensity fluctuations, 
we need angular resolution ^res corresponding to Rs- In the Einstein-de Sitter universe, this 
angular resolution is given by 

i?s R. ^26) 



2(yrT^-i) ■ Vio-3 

where the last value is for 1 + z = 10. 



^ arcmin, (27) 



4. Intensity fluctuations due to nonlinear clumps 

If we observe with high enough angular resolution, we should be able to detect features from 
individual clumps, rather than from linear matter inhomogeneities. More concentrated objects 
will be potentially stronger sources. 

Hierarchical clustering model of structure formation in the universe tells us how nonlinear 
clumps formed in the process of gravitational instability growing. Let us estimate typical values 
of the comoving size of clumps -Rtyp and the mean separation among clumps -Rtyp at 1 + z = 10. 

The comoving number density of clumps with in the mass range M ~ M + dM at z is given 
by the Press-Schechter (1974) formalism: 



N{M, z)dM 



2 po <Jc 



vr M D{z) 



1 da{M) 



a'^{M) dM 



exp 



2a^{M)D^{z) 



dM, 



(28) 
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where pQ is the mean comoving density of matter, 5c is the overdensity threshold, D{z) is the 
growing function of perturbations, and a{M) is the variance of density field smoothed with mass 
scale M. Since a{M) is a decreasing function of M, the comoving number density (equation 
exponentially falls down at large mass. Therefore we can estimate the typical mass of clumps Mtyp 
by a relation: 

Since the density of virialized objects is ISvr^ times the mean density of the universe in spherical 
collapse model, i?typ is related to Mtyp by i?typ = (3Mtyp/47r/9o)"^''^(187r^)~^/^. Furthermore, the 
mean comoving separation Rtyp = n~^/^ is given from the comoving number density: 

POO 

n{Rtyp,z)= N{M,z)dM. (30) 

We make use of the fitting formula for (t{M) in CDM models by ( [Kitayama and Suto 199^ ): 

o-(M) oc (1 + 2.208mP - 0.7668m2p + 0.7949m3P)- i , (31) 

which gives a good fit between lO"'^ < m < 10^. Here p = 0.0873, and m = M(Oo/i^) 7(10^^^©). 
In the standard CDM model (JIq = 1) with h = 0.5 and fis = 1, the above equations yield 

i?typ = 0.015/1"^ Mpc, R^ = 0.4/1-1 Mpc, (32) 

at 1 + z = 10. In terms of Az^/z/, these numbers correspond to 2 x 10^^ and 4 x 10"^. 
The one-dimensional velocity dispersion of these typical clumps is estimated as 



i;typ = VG'Aftyp(l + z)/(2i?typ) = 30kms-\ (33) 
which corresponds to IS.v/v = 1 x 10^^. Hence the effective overdensity of the clump is 

Jeff = ISvr^ X ^(iJtyp)/!^ = 30. (34) 

In CI (609 ^m) case, for example, bumps of A///cmb = 1 x 10"^ are expected to appear with mean 
interval of 4 x 10~^ in /S.v/v, if observation is made with frequency resolution /S.v/v = 1 x 10^"^ 
and angular resolution O^^s — 1 arcsec. Here the angular resolution is related through equation ( p6| ) 
with the clump size Rtyp- In more conventional units in radio astronomy, this line has intensity 
of ATiinc'Utyp = 0.008 K km s~i, where ATiinc is the increment of brightness temperature. In CI 
(370 /^m) case, the amplitude of intensity fluctuations are estimated to be A///cmb = 4 x 10"*^, 
that is ATiine'Wtyp = 0.02Kkms~i, with the same resolution. Alternatively, if we assume that most 
carbon is in the form of CII, the signal of CII (158 /im) line is estimated to be A///cmb = 3 x 10~^, 
that is ATiinet'typ = 0.07Kkms~i. 
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5. Characteristics of individual lines 

The observed signal will depend upon the ionization state and abundance of each element. 
Here, we present a simple consideration regarding the ionization state based on ionization 
potential. More realistic treatment of ionization of each element requires detailed computation 
including temperature and density inhomogeneities of gas and radiative transfer, which is beyond 
the scope of this paper. 

CI 609 //m and 370 /im lines have the longest wavelength among the lines discussed here. 
If 1 + ~ 10-20, they are easiest to observe from ground. In addition, they can easily be 
distinguished from others. Even if the signal is weak, one can extract the signal by cross-correlating 
different portions of observed spectrum corresponding to these two wavelengths. The most serious 
concern is that the same star that produce carbon, will quickly ionize much of it to CII. The 
ionization potential of CI, 11.3 eV, is smaller than that of hydrogen, so we may have large mostly 
neutral regions containing CII. 

Since carbon is so easily ionized, CII 158 /im line may be a more promising candidate 
to be detected. CII is the brightest submillimeter line in our galaxy ( Mizutani et al. 1994| ). 



Unfortunately, this is a single line and its higher frequency makes ground-based detection more 
challenging. 

Other promising candidates are Nil 204 /im and 122 //m, and 01 146 fim lines. Nitrogen 
is expected to be abundant only after a relatively low redshift, since it is mainly produced in 
metal-enriched stars through incomplete CNO burning. Moreover, the ionization potential of 
nitrogen is 14.5 eV, so it is more difficult to ionize than hydrogen. On the other hand, 01 
abundance is supposed to peak around the reionization epoch; the ionization potential of oxygen 
is 13.6 eV, nearly equal to that of hydrogen, and charge exchange reaction (O^ + H <-> O + H"*") 
imply that noi/noii ~ nm/nmi- 

The other, relatively short-wavelength lines are not likely to be observed from the ground, 
unless there was significant amount of metal production by 1 + z ~ 25-50. Moreover, NIII and 
OIII probably never have been dominant species in each element; the ionization potential of Nil 
and on are 29.6 eV and 35.1 eV, respectively. However, their line intensity is so strong that 
they may be useful to place an interesting constraint on the state of matter and radiation field at 
1 + z ~ 25-50. 



6. Alternative estimate of CII intensity 

We can make an alternative estimate of the CII 158/xm line intensity, clued by the 
measurement that roughly 1 percent of the luminosity from low redshift galaxies is emitted in 
the CII line ([Malhotra et al. 1997] ). Though there is a possibility that the metal contamination 
may have somewhat changed the fraction, we assume that the same fraction (1 percent) of the 
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luminosity from small cloudlets at z = Zg is emitted in the CII line. 

We estimate the luminosity density £cii of the CII emission instead of /ipt'iine'^M^Mi in 
equation (^). We assume that the small cloudlets evolve with a constant star formation rate B 
per comoving unit volume from the Salpeter type initial mass function (j){M) (Salpeter 1955) and 



com 



achieve metallicity Z{zs) = O.OIZq by z = Zg. We use the integrated abundance yields of type II 
supernovae tabulated by Thielemann, Nomoto, and Hashimoto (1993) to derive the star formation 
rate 

° ^ V 10 / V0.5y V2.1 X 10-7cm-3y V 0.0043 J \ 0.01 J ' ^ ^ 

where n^™ is the comoving number density of hydrogen. Then the comoving luminosity density 
in the B-band is given by 

= r^"s=°"^t„(zs)(/>(M)Lij(M)dM+ / B''°'^T{M)(f>{M)LBiM)dM, (36) 

J JMr 

where tu{zg) is the age of the universe at z = Zs, t{M) and Lb{M) are the age and the 
B-band luminosity of a star whose mass is M, respectively, and M^- is the mass which satisfies 

T{Mr) = tu{Zs). 

The resultant comoving luminosity density is >C™jf = O.Ol/Zg'™ = 4 x lO''L0Mpc~^ (at 
1 + 2:3 = 10), 7 X IO'^LqMpc"^ (at 1 + = 20). We find that these are in good agreement with 
the previous estimate h-pv\inen'^^ ^ui = 6 x lO''L0Mpc~'^. We conclude that the CII intensity 
fluctuations estimated alternatively here are almost the same magnitude as the previous result 
(table I). 



7. Conclusion 

We have proposed a new and promising method to detect the first objects: search for carbon, 
nitrogen and oxygen line emission at 1 + z = 10-20. We conclude that the intensity fluctuations 
originated from the emission are detectable in radio band (0.5 mm-1.2 cm). If observation is carried 
out with angular resolution ^res = 1 arcmin and spectral resolution Az^/i/ = 10~^, the magnitude 
of the smoothed intensity fluctuations corresponding with linear matter density fluctuations is 
predicted as Alnne/Ztot = (1 and 3) x 10^® (in CI 609 and 370 //m cases) from sources at 1 + z = 10, 
in the standard CDM model. The predictions for other lines are A7iino/-^tot = (20, 70, 3, 20) x 10^^, 
respectively, in [158(CII), 146(01), 204(NII), 122(NII)] /im cases. 

Observation with ^rcs = larcsec and /S.u/u = 10~^ resolve individual protogalaxies with 
velocity dispersion of 30kms~^. We predict line intensities relative to the microwave background of 
A/iine/-ftot = (1, 4, 30, 80, 4, 30) x lO""^ or equivalently ATuneWtyp = (0.8, 2, 7, 20, 1, 6) x 10"^ K kms"^ 
for [609(CI), 370(CI), 158(CII), 146(01), 204(NII), 122(NII)] /im lines, respectively. If detected. 
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these features will provide useful information on matter inhomogeneities at 1 + z = 10-20. They 
will also probe chemical evolution and reionization time Zj-c itself. 

We thank John H. Black, Bruce T. Draine, Thomas Herbig, Ryoichi Nishi, Jeremiah P. 
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Fig. 1. — Smoothed intensity fluctuations: Solid and Dotted lines show the signals for CI (609 fim) 
and CI (370 /um) emission, respectively, if most of the carbon is in the form of CI. Broken line shows 
the signal for CII (158 jum) emission in the case where CII ion state is preferred. 



-14- 



Peebles, P. J. E. 1993, Principles of Physical Cosmology: Princeton Univ. Press 
Press, W. H., and Schechter, P. 1974, ApJ, 187, 425 
Salpeter, E., E. 1955, ApJ, 121, 161 

Schneider, D. P., Schmidt, M., and Gunn, J. E. 1991, AJ, 102, 837 

Steidel, C. C, Giavalisco, M., Pettini, M., Dickinson, M., and Adelberger, K. L., 1996, ApJ, 462, 
L17 

Stutzki, J., ct al. 1997, ApJ, 477, L33 

Tegmark, M., Silk, J., Rees, M. J., Blanchard, A., Abel, T., and Palla, F. 1997, ApJ, 474, 1 

Thielemann, F.-K., Nomoto, K., and Hashimoto M. 1993, Origin and Evolution of the Elements, 
N. Prantzos, E. Vangioni-flam and M. Casse: Cambridge, 297 

Woosley, S. E., and Weaver, T. A. 1995, ApJS, 101, 181 

Yamamoto, S., and Saito, S. 1991, ApJ, 370, L103 

Young, K., Cox, P., Huggins, P. J., Forveille, T., and Bachiller, R. 1997, ApJ, 482, LlOl 



This preprint was prepared with the AAS lAT^^X macros v4.0. 



-15- 



Table 1: Transition lines and intensity fluctuations: the estimates of A/uneZ-ftot assume that the 
each ion is the preferred ionization state. 



Species 


Transition 


Wavelength 




A-^lineZ-^tot 






(urn) 


(s-i) 




CO 


'Pi ^ 'Po 


609 


7.93 X 10-s 


1 X 10-6 


CO 


^ 'Pi 


370 


2.68 X 10-"^ 


3 X 10^6 


C+ 


2 P _^ 2 p 
^3/2 ~^ ^1/2 


158 


2.36 X 10-6 


2 X 10-5 


N+ 


'Pi - 'Po 


204 


2.13 X 10-6 


3 X 10-6 


N+ 


'P2 'Pi 


122 


7.48 X 10-6 


2 X 10-5 


N++ 


2 p , ^ 2 p , 
^3/2 ^1/2 


57.3 


4.77 X 10"^ 


2 X 10-''^ 


0° 


'Pi - 'P2 


63.2 


8.95 X 10"^ 


9 X IQ-^ 


0° 


'Po - -^^1 


146 


1.70 X 10-^ 


7 X 10-5 


0++ 


'P2 ^ 'Pi 


51.8 


9.75 X 10-^ 


5 X 10-2 


0++ 


'Pi ^ 'Po 


88.4 


2.62 X 10-5 


7 X 10-^ 



